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Abstract 

Direct CP violation in the hadronic decays is investigated near 

the peak of the taking into account the effect of p—LV mixing. Branching ratios for 
processes 5±,o 

—> and B —> ujtt are calculated as well. We find that the 

CP violating asymmetry is strongly dependent on the CKM matrix elements. For 
a fixed , the CP violating asymmetry, acp, has a maximum of order —40% 
to —70% for B^ —> p^{u;)Tr^ when the invariant mass of the tt+tt” pair is in the 
vicinity of the u resonance. The sensitivity of the asymmetry to is small in 

that case. Moreover, we find that in the range of which is allowed by the 

most recent experimental branching ratios from the BABAR, BELLE and CLEO 
Collaborations, the sign of sin (5 is always positive. Thus, a measurement of direct 
CP violation in decays B^ —> Tr^Tr^vr^ would remove the mod(7r) ambiguity in the 
determination of the CP violating phase angle a. 
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1 Introduction 


In the Standard Model, CP violating phenomena arise from a non-zero weak phase angle 
in a complex matrix allowing flavour violation in the weak interaction: the Cabbibo- 
Kobayashi-Maskawa (CKM) matrix. Although the source of CP violation has not been 
well understood up to now, physicists are striving to increase their knowledge of the 
mechanism. Many theoretical studies |2] (within and beyond the Standard Model) 
and experimental investigations have been conducted since the discovery oi CP violation 
in neutral Kaon decays (1964). According to theoretical predictions, large CP violating 
effects may be expected in B meson decays. In the past few years, several facilities 
have started to collect events on B decays and most of them refer to branching ratios. 
Generally, the main theoretical uncertainties apart from the CKM matrix elements are the 
hadronic matrix elements, where non-factorizable effects are involved. As regards hadronic 
matrix elements and non-factorizable effects, a new QCD factorization approach [H] has 
been proposed. This QCD factorization approach includes all radiative diagrams (gluon 
exchange) but will not be the subject of this paper. For the CKM matrix elements, 
uncertainties in the parameters p and rj have been reduced and this allows us to predict 
CP violating asymmetry in B decays more accurately than before. This will give us an 
excellent test for the Standard Model and may lead to suggestions of new physics. 

Direct CP violating asymmetries in B decays occur through the interference of at least 
two amplitudes with different weak phase (j) and strong phase 6. In order to extract 
the weak phase (which is determined by the CKM matrix elements), one must know 
the strong phase 6 and this is usually not well determined. In addition, in order to 
have a large signal, we have to appeal to some phenomenological mechanism to obtain 
a large 6. The charge symmetry violating mixing between and a; can be extremely 
important in this regard. In particular, it can lead to a large CP violation in B decay 
such as B^ —p°(a;)7r° —>• 7r+7r“7r°, because the strong phase passes through 90" at the uj 
resonance m 0 El 

We have collected all the recent data for b to d transitions, but we shall focus on the 
CLEO, BABAR and BELLE branching ratio results. We also shall use the latest values 
for CKM parameters. A, A, p, and p. The aim of the present work is to constrain the CP 
violating calculation in B^ —> p^{uj)ti^ —>■ 7r’''7r“7r‘^, including p — uj mixing and using the 
most recent experimental data for the branching ratios for B ^ pn decays. In order to 
extract the strong phase 6, we use the naive factorization approach, in which the hadronic 
matrix elements of operators are saturated by vacuum intermediate states. Moreover, we 
approximate non-factorizable effects by introducing an effective number of colours, N^f, 
In this paper, we investigate hve phenomenological models with different weak form factors 
and determine the CP violating asymmetry for B^ p^{uj)'K^ —>■ in these models. 

We select models which are consistent with all the latest data and determine the allowed 
range for N^f (1.09(1.11) < N^f < 1.68(1.80)). Then, we study the sign of sin 5 in this 
range of N^f in all these models. We also discuss the model dependence of our results 
in detail. 

This paper is structured as follows. In Section 2, we introduce the effective Hamiltonian 
based on the Operator Product Expansion (OPE) including Wilson coefficients. We also 
present the formalism of p — a; mixing and its application to the CP violating asymmetry 
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in decay processes. In Section 3, the CKM matrix and the relevant form factors are 
discussed. In Section 4, we present numerical results for the CP violating asymmetry in 
—> 7r^7r“7r‘’ which is followed by discussion of these results. In Section 5, branching 
ratios for decays such as —>■ and B~ —>■ a;7r“ are investigated. From the 

CLEO, BABAR and BELLE experimental data for these branching ratios, we extract the 
range of N^f allowed in these processes and the results are also discussed. In the hnal 
section, we summarize our results. Comments on form factors, CKM matrix parameter 
values, p, p, and conclusions are also given in this section. 


2 CP violation in —> tt+tt tt® 


2.1 Effective theory 

In any phenomenological treatment of the weak decays of hadrons, the starting point is 
the weak effective Hamiltonian at low energy [ 7 j from which, the decay amplitude can be 
expressed as follows. 


A{B PV) = 


% v„,v:^(c,{PV\o';\B) + C2 {PV\ 0'i\B)) 

10 

-Vu,v:,Y.c,{pv\oab) 

i=3 


+ h.c. , (1) 


where {PV\Oi\B) are the hadronic matrix elements. They describe the transition between 
initial and hnal states with the operator renormalized at scale /i and include, up to now, 
the main uncertainties in the calculation since they involve non-perturbative effects. Gp is 
the Fermi constant, Vqkm is the CKM matrix element, Ci{fi) are the Wilson coefficients, 
Oi{fi) are the operators from OPE jH]. The operators Oi, the local operators which govern 
weak decays can be written as. 


0“ = - j5)Ul3U0-f^{l - 75)&a , 

O3 = 97^(1 - 75)^ ~ , 

gf 

O 5 = g7M(i - 75 )^ , 

q' 

Or = ^g 7 M(l -75)^5^ e 5 /gV(l + 75)9' , 

q' 

O 9 = -75)&X]e5/gV(l -75)g , 


O2 = ?7m(1 - 75)mm7^(1 - 75 )^ , 

O 4 = - 75)&p - 75)ga , 


q' 

0 & = qal^,{^ - 75 )&p X] 9 p 7 ''(l + lb)q'a > 

q' 

3 

Os = -ga7M(l - 75)&P + l3)q'a , 

q' 

3 

Oio = -ga7M(l - 75)&P V^p7''(l - 75)g« , 

q' 


( 2 ) 


where q' = u, d, s, c and e,/ denotes its electric charge. As regards the Wilson coeffi¬ 
cients Eiiiniiniini, they represent the physical contributions from scales higher than p. 
Since QCD has the property of asymptotic freedom, they can be calculated in perturba¬ 
tion theory. Usually, the scale p is chosen to be of order 0{mb) for B decays and Wilson 
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coefficients have been calculated to the next-to-leading order (NLO). For more details see 

Ref. |T!^ . 


2.2 p — uj Mixing 

Let A be the amplitude for the decay B —>■ p^7i —> then one has, 

A = -|- (7r7r“7r'''|iL‘^|i?) , (3) 

with and being the Hamiltonians for the tree and penguin operators. We can 
dehne the relative magnitude and phases between these two contributions as follows, 

A = (7r7r-7r+|i/^|R)[l + re*V‘^] , 

A = (7r7r+7r-jH^jB} [1 + , (4) 


where 6 and 0 are the strong and weak phases, respectively. The phase 0 arises from the 
appropriate combination of CKM matrix elements, and, assuming top quark dominance, 
0 = arg[(Vi6V^^)/(14bV^rf)]- As a result, sin0 is equal to sin a, with a dehned in the 
standard way [2]. The parameter, r, is the absolute value of the ratio of tree and 
penguin amplitudes: 

(pOTrjll^jB} ■ ^ ^ 

In order to obtain a large signal for direct CP violation, we need some mechanism to 
make both sin 6 and r large. We stress that p — uj mixing usi has the dual advantages 
that the strong phase difference is large (passing through 90" at the uj resonance) and 
well known isi ini- With this mechanism, to hrst order in isospin violation, we have the 
following results when the invariant mass of 7r’''7r“ is near the uj resonance mass, 

(7r7r"7r+|iL^|R) = + —tp , 

SpSco Sp 

{n7r~n^\H^\B) = -^^Upu^Puj + —Pp ■ ( 6 ) 

SpSuj Sp 

Here tyiV = p or a;) is the tree amplitude and pv the penguin amplitude for producing 
a vector meson, V, Pp is the coupling for tt+tt”, Hp^ is the effective p — uj mixing 

amplitude, and sy is from the inverse propagator of the vector meson V, 

Sy = s - rriy + imyVy , (7) 


with a/s being the invariant mass of the tt+tt” pair. We stress that the direct coupling 
UJ TT+TT” is effectively absorbed into Hp^ [2], leading to the explicit s dependence of 
Hp^. Making the expansion np^(s) = np^(m^) + (s - ''7i^)np^(m^), the p - a; mixing 
parameters were determined in the £t of Gardner and O’Connell |n|: 3fJe np^(m^) = 
-3500 ± 300 MeV^ Am tip^{ml) = -300 ± 300 MeV^ and ti'p^iml) = 0.03 ± 0.04. In 
practice, the effect of the derivative term is negligible. From Eqs. (0110 one has, 


re 


npcjPa; “1“ ^ujVp 


( 8 ) 


Dehning, 
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(9) 


— = , — = ae“" , 



where 6a, 6^, and 6q are strong phases (absorptive part). Snbstitnting Eq. Q into Eq. (jH)), 
one hnds, 






Suj + 




iSoc 


( 10 ) 


0 ( 6 *'^“, , and will be calcnlated later. In order to get the CP violating asymmetry 

acp, sin0 and cos0 are needed, where 0 is determined by the CKM matrix elements. In 
the Wolfenstein parametrization, the weak phase comes from [VtbV*^/VubV*^] and one has 
for the decay B —>• p(a;)7r, 


sin (p 


COS0 


_ V _ 

V[p(i -p) -v‘^? + v^ ’ 

p(l - p)-r[^ 

V'[p(l -p) + 


( 11 ) 


The valnes used for p and p will be discussed in Section Id. 11 With the decay amplitude 
given in Eq. dU, we are ready to evaluate the matrix elements for In the 

factorization approximation CHI, either p°(a;) or is generated by one current which 
has the appropriate quantum numbers in the Hamiltonian. Eor these decay processes, 
two kinds of matrix element products are involved after factorization; schematically (i.e. 
omitting Dirac matrices and colour labels) one has (p°(a;)|('UM)|0)(7r=*"P|(-u6)|i?^P) and 
( 7 r^P|(gig 2 )| 0 )(p°(a;)I(126)11?^’°) with qi{i = 1,2) = u,d. We will calculate them in some 
phenomenological quark models. 

The matrix elements for B ^ X and B —>■ X* (where X and X* denote pseudoscalar and 
vector mesons, respectively) can be decomposed as follows [T?l] . 


{X\J,\B) = [pb+Px 


m 


B 


m 






k] Fi(p2) + 




( 12 ) 


and 




mB + mx* 
e* ■ k 


ms + mx* 


(^^upae*''pBPx*^{k ) + *<( Cimp + mx*)Ai{}i^) 


{Pb + Px*)p,A2{}A) -p—2mx* • k^Aj,{lC‘) 


+ i^-^2mx* ■ kf,Ao{k^) , (13) 


where Jp(= g 7 ^(l — 75 )^) is the weak current with q = u,d, k = Pb — Px(x*), and is 
the polarization vector of X*. Fq and Fi are the form factors related to the transition 
0“ —0“ and Hq, Hi, 2 I 2 , 2 I 3 , and V are the form factors which describe the transition 


4 













0 —1 . Finally, in order to cancel the poles at = 0, the form factors must respect 

the constraints: 

Fi(0)=Fo(0), A,{0)=Ao{0). (14) 

They also satisfy the following relations: 


Asik^) 


ms + mx* 
2mx* 


Ai{e) 


ms — mx* 
2mx* 


A2{k^) . 


(15) 


By using the decomposition in Eqs. (HI d, one obtains the following tree operator 
contribution for the process > p‘^(a;)7r‘^: 


tp = mB\Pp 


{C[ + 


— 77 C 2 ) 
N-ff 2 ^ 


fpFi{ml) + UAo{ml) 


(16) 


where fp and fn are the decay constants of p and tt, respectively, and C[ are the Wilson 
coefficients with values listed in Table [H We hnd ^ tp, so that 




-fpFiiml) + UAo{ml) 


fpFi{ml) + f^Ao{ml) 
After calculating the penguin operator contributions, one has, 

Pw 






VaV,", 




(17) 


(18) 


and 


_ b\Pp\ I _ C'3)[/pFi(m^) + /^Ao(m^)] 

Poj I iVc 

- \m + ^C') - (Q + ArrCmUMK 




+ 2^^^^ ^e/Z^s) + (^9 + ^eff ^lo)]/p-^l(’^p) 

\(r' J. _ \(r' ^ [ 2mlUAo{ml) 

® ® 2 ® ^ [{md +md){mb +rrid) 

+ 2^^10 + UAo{ml)] )■ , (19) 


where Pp is the c.m. momentum of the decay process. In Eqs. (IT^ I19|l . is written as. 
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Plo = rnslPp 


(Q++ (cj+ 


1 


fpFiiml) 


^^7 + + (C'g + 


/p^l(7«J) 


2 ^^8 + - (C'e + 


2mlf^Ao{ml) 

{md + md){mb + rud) 

- (^4 + -^C',) [UAo{ml) + /pFi(mJ)] 


mff 


f„Ao(ml) + 7;fpFi{'mp) 


, ( 20 ) 


and the CKM amplitude entering the b ^ d transition is, 


VtbV,:, 

\l{^-p)^ + v^ r A^y^ 

sin 7 

Vu,v:, 

(l-AV2)Vp2 + p2 V 

sin f3 


with (3 and 7 dehned in the unitarity triangle as usual. 


( 21 ) 


3 Numerical inputs 

3.1 CKM values and quark masses 

In our numerical calculations we have several parameters: and the CKM matrix 

elements in the Wolfenstein parametrization. The CKM matrix, which should be de¬ 
termined from experimental data, is expressed in terms of the Wolfenstein parameters, 
A, A, p, and 7 [20]. Here we shall use the latest values [23 which have been extracted 
from charmless semileptonic B decays (|I4fe|), charmed semileptonic B decays (|Wb|), s 
and d mass oscillations and CP violation in the kaon system (p, 7 ): 

A = 0.2237 , A = 0.8113 , 0.190 < p < 0.268 , 0.284 < p < 0.366 . (22) 

These values respect the unitarity triangle as well. The running quark masses are used 
in order to calculate the matrix elements of penguin operators. The quark mass is taken 
at the scale p ~ m;, in H decays. Therefore one has [221, 

mu(p = mb) = 2.3 MeV , mb{p = mb) = 4.9 GeV , md{p = mb) = 4.6 MeV , (23) 

which corresponds to ms{fi = 1 GeV) = 140 MeV. As regards meson masses, we shall use 
the following values HH: 

mso = 5.279 GeV , m^± = 0.139 GeV , m^^o = 0.135 GeV , 

mpo = 0.769 GeV , m^ = 0.782 GeV . (24) 
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3.2 Form factors and decay constants 

The form factors Fi{k‘^) and Aj{k‘^) depend on the inner strnctnre of hadrons. In order 
to gange the model dependence of the resnlts, we will adopt three different theoretical 
approaches. The first was proposed by Baner, Stech, and Wirbel jlH] (BSW model). They 
nsed the overlap integrals of wave fnnctions in order to evalnate the meson-meson matrix 
elements of the corresponding cnrrent. The second approach was developed by Gno and 
Hnang (GH model) [2H]- They modihed the BSW model by nsing some wave fnnctions 
described in the light-cone framework. The last model was given by Ball 121 ESI- In this 
case, the form factors are calcnlated from QGD snm rules on the light-cone and leading 
twist contributions, radiative corrections and S't/(3)-breaking effects are included. The 
explicit k"^ dependence of the form factors is uniEi, 


and [21 121 El 


Fi(fc2) = 


hi 


m? 


Ao{k^) = 






hi 


l-d,A^ + bi (K) 

\^bJ 


2 ’ 


Ao{k^) = 


h 


Ao 




2 ’ 


(25) 


where n = 1,2, and and mi are the pole masses associated with the transition current, 
hi and h^^ are the values of the corresponding form factors at = 0, and di, bi (z = 0,1) 
are parameters in the model of Ball. In Table El we list the relevant form factor values 
at zero momentum transfer m El El El El for B ^ 71 and B ^ p transitions. The 
different models are defined as follows: models (1) and (3) are the BSW models where 
the dependence of the form factors is described by a single and a double-pole ansatz, 
respectively. Models (2) and (4) are the GH model with the same momentum dependence 
as models (1) and (3). Finally, model (5) refers to the Ball model. We define the decay 
constants for pseudo-scalar (/p) and vector (/y) mesons as usual by. 


(^(?) 1 ^ 17 / 275 ^ 210 ) = -ifpqi, , 

v^(F(g) I hi 7/2^210) = fvmvev , (26) 

with being the momentum of the pseudo-scalar meson, and my and ey being the mass 
and polarization vector of the vector meson, respectively. In our calculations we take m 

U = 132 MeV , fp-U = 221 MeV . (27) 

In practise the p and uj decay constants are very close, and as a simplihcation (with little 
effect on the results), we chose fp = fu}- 


4 Results and discussions 

A previous analysis j2H] has been conducted showing the dependence on the GKM matrix 
elements and form factors of the direct CP violating asymmetry. Here, we update our 
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investigation by taking into account the latest values of the Wolfenstein CKM parameters, 
p and rj, and also by analysing more B decays. In the following numerical calculations, we 
apply the formalism detailed previously and investigate more precisely. We 

hnd that for a hxed there is a maximum value, amax, for the CP violating parameter, 
acp, when the invariant mass of the 7r’''7r“ pair is in the vicinity of the ut resonance. 
In Figs. ^ and 121 CP violating asymmetries for B^ —> for k'^/ml = 0.3 with 

NCf = 1.09(1.68), and k'^/ml = 0.5 with N/P = 1.11(1.80), are plotted, respectively, 
and for limiting values of CKM matrix elements. Graphic results are shown only for the 
model (1), as an example. We have investigated hve models, with hve different form 
factors in order to test the model dependence of acp- 

Concerning the maximum CP violating asymmetry for B^ amaxi h varies 

from —51%(—38%) to —84%(—69%) in the allowed range of p,T] for k'^/ml = 0.3(0.5). 
From the numerical results listed in Table El for = 1.09(1.11) and = 

1.68(1.80), we can see that the hve models fall into two classes: models (1,3) and (5) 
and models (2) and (4). For models (1,3) and (5), and for = 1.09(1.11), the 

maximum asymmetry, amax-, is around —54%(—40%) for the set {pmax,Vmax) and around 
—69%(—53.6%) for the set {pmimVmin)-, leading to the ratio between them being around 
1.28(1.34). In each of these models and for = 1.68(1.80), the maximum value of 

the asymmetry, amax, varies from —62.6%(—48.6%) for the set {pmax,Vmax) to around 
—77.3%(—64.6%) for the set {pmin,Vmin)- In that case, the ratio is equal to 1.23(1.32). If 
we consider models (2) and (4), the maximum asymmetry, amax, where = 1.09(1.11), 
is around —63.5%(—48%) for the set {pmax,Vmax) and around —78.5%(—62%) for the set 
{Pmin,riniin)- This yields aratio 1.24(1.29). When = 1.68(1.80), one has a maximum 

asymmetry around —71%(—56.5%) for the set {pmax,Vmax) and around —84%(—69%) for 
the set {pmin,Vmin), leading to a ratio around 1.18(1.22). 

From all these results, many comments can be enumerated. Although the maximum 
asymmetry, amax, still varies over some range in the B^ —> 7r+7r“7r° decay, we stress that 
by using more accurate CKM element values than before, a more precise CP violating 
asymmetry is obtained. The reason is primarily the matrix elements Vtd and Wfe which are 
involved in the b ^ d transition through the ratio of to tp. In our previous CP violation 
study I2S] for the process B —>■ tt , we found that the ratio between the maximum 
and minimum asymmetry, related to the minimum and maximum set of (p, p), was around 
1.6. By comparison, in the present work, this ratio is reduced to 1.3. The difference is 
related to the improvement in the measurement of the CKM matrix elements, and shows 
the strong effect of the CKM parameters, p and p, on limiting asymmetry values. 

With regard to the CKM matrix elements, it appears that if we take their upper limit, 
we obtain a smaller asymmetry, acp, and vice-versa. As we found before, there is still a 
strong dependence of the CP violating asymmetry on the form factors. The difference 
between the two classes of models, (1, 3, 5) and (2,4), comes mainly from the magnitudes 
of the form factors. In fact, the form factor Fi{k‘^), which describes the transition B —>■ 7i, 
is mainly responsible for this dependence. In both classes, we hnd a stronger dependence 
of the CP violating asymmetry on the CKM matrix elements than that on the form 
factors or the effective parameter N/P. The difference observed in our results between 
k"^/ml = 0.3 and k"^/ml = 0.5 arises from the k‘^ dependence of the Wilson coefficients in 
the weak effective Hamiltonian. Finally, since N/P (treated as a free parameter) is related 


to hadronization effects through the factorization approach, it is not possible to determine 
its value accurately (since non-factorizable effects are not well known). That is why the 
asymmetry also varies in some range of . It is obvious that a more accurate value for 
(which requires a more accurate approach with non-factorizable effects being taken 
into account), and hadronic decay form factors (which requires better understanding for 
pionic structure and the i? —> tt transition) are needed in order to determine the CKM 
matrix elements. 

In spite of all the uncertainties mentioned above, we stress that the p — uj mixing 
mechanism in the B ^ pn decay can be used to remove ambiguity concerning the sign 
of sin 5. As the internal top quark dominates the 6 —>■ d transition, the weak phase 
in the asymmetry is proportional to sin a (= sin0), where a = arg . Hence 

ud III) _ 

knowing the sign of sin 6 enables us to determine that of sin a from a measurement of the 
asymmetry, acp- In Fig. El we show sinh as a function of for —>■ when 

we have maximum CP violation. Then, in our determined range of N^f, (1.09(1.11) < 
< 1.68(1.80)), one hnds that its sign is always positive for all the models studied 
and for all the form factors. Therefore, by measuring the CP violating asymmetry in 
B^ —>■ 7T~^7r~7T^, we can remove the mod(7r) ambiguity which appears in the determination 
for a from the usual indirect measurements which yield sin2Q!. In Fig. EJ the ratio of the 
penguin and tree amplitudes, as a function of , is plotted for limiting values of the 
CKM matrix elements, p, p, for the process B^ —>■ n~^7r~7T^. Even though one gets a larger 
value of sind around = 1, for B^ —>■ without p — uj mixing, one still has a 

small value for r around this value of . In that case, the CP violating asymmetry, 
acp, remains very small without p — uj mixing. 


5 Branching ratios for 

5.1 Formalism 

The direct B —>■ transition is the main contribution to the decay rate. In our case, 
to be consistent, we should also take into account the p — uj mixing contribution to the 
branching ratio, since we are working to the hrst order of isospin violation. The derivation 
is straightforward and we obtain the following form for the branching ratio for B p^n: 


BR{B p°7r) 



) Oio) 

+ 

02 ) - , aio) 

^^pLO 

{Sp - ml) + imj:^ 


(28) 


In Eq. (j28|l Cp is the Fermi constant, F^ is the B total decay width, and is an 
integer related to the given decay. Ay and Ay are the tree and penguin amplitudes, and 
I//’, Vf represent the CKM matrix elements involved in the tree and penguin diagrams, 
respectively: 

Vl=\VM\, and Vf = \VM\. (29) 
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The effective parameters, a*, which are involved in the decay amplitude, are the following 
combinations of effective Wilson coefficients: 


a2j — C2j + 




^C' 


a 2 j-i = CL_i + for j = 1, • • • , 5 . (30) 


n: 


5.2 Calculational details 

In this section, we give full details of the theoretical decay amplitudes for decays involving 
the b to d transition. Two of these decays involve p — uj mixing. They are B~ p^7T~ 
and —>■ p° 7 r°. The other two decays are B^ —>■ p~Ti^ and B~ p~tt^. We list in the 
following, the tree and penguin amplitudes which appear in the given transitions. 

For the decay B~ = 32 in Eq. 


v^/ 1 T(oi, 02 ) = + a2/..4o(mJ) , 


(31) 


3 1 

V2Ap{a3i • ■ ■ , Oio) = fpFi(m‘^p) —04 + 2 + ^ 9 ) + 2 ®^° 


+ /7r^o(^^)s ~ 2 (a 6 + Og) 


mz 


_{mu + md){mb + mu) 


+ Oio f ; (32) 


for the decay B —ujir (a^ = 32 in Eq. (j2Hl)), 
V2Al{ai,a2) = ai/pFi(mJ) + a2UAo{ml) , 


(33) 


■ ■ ■ )®io) — fpFi{m'l)l 2(03 + Os) + -(a? + 09 ) + (^4 — x^io) 


+ fTrAo{ml.)< —2{ae + Og) 


m: 


_{mu + md){mb + m„)_ 


+ 04 + Oio > ; (34) 


for the decay B^ —> p°7r° {ak = 64 in Eq. (I^ i. 

2^p ( 01 , 02 ) = Oi/pFi(mJ) + aiUAo{ml) , (35) 


2 ^p ( 03 , • • • , Oio) — fpFi{mp)'^ —04 + -( 3 a 7 + 3 a 9 + Oio) [■ + 


/jr74o(m^)< —04 + (200 — Og) 


mZ 


_2md{mb + md) 
for the decay B^ {ak = 64 in Eq. (I^ i. 


+ + Sog + oio) ^ ; (36) 


2A^(ai,a2) = -aifpFi{m‘^p)+aif^Ao{ml) , 


(37) 
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2^0) (®3) ■ ■ ■ ) ®lo) — fpPl{^n) 1 ~2(a3 + Os) — 04 — -(Oy + Og — Oig) 


+ fnAo{m‘^) < —04 + (206 — Og) 


m: 


_2md{mb + md) 

for the decay —>• p“7r+ (a^ = 16 in Eq. (I^ i. 


+ + 3®9 + ®io) f ; (38) 


A J ( 01 , 02 ) = a 2 fpFi{ml) 


(39) 


(® 3 ) • • •) Olo) — (04 + oio)/p-f"i(’^p) ; 
for the decay B~ —>■ p~ (a*. = 32 in Eq. (OHjl i. 


v/2AJ’(oi, 02) = 02/pEi(m^) + oi/^Ao(m^) , 
■\/2^p (03, • • • , Oio) = /pEi(m^)(o4 + Oio) + 


UAo{ml 


-“4 - 2 ( 3^7 


SOg — Oio) + (2 o 6 — Og) 


mf 


_2md{mb + mrf)_ 


(40) 


(41) 


(42) 


Moreover, we can calculate the ratio between two branching ratios, namely BR{B^ —> 
p^TT^) and BR{B^ —>■ p^vr^), in which the uncertainty caused by many systematic errors 
is removed. We dehne the ratio, R^,, as: 


BR{B^ p±7r^) 
BR{B^ p07r±) ■ 


(43) 


5.3 Numerical results 

The numerical values for the CKM matrix elements V^' , p — oj mixing amplitude 
and particle masses mv,p, which appear in Eq. (I28j) . have been reported in Sections 12.21 
and El The Fermi constant is taken to be Gp = 1.166391 x 10“^GeV“^ [H], and for the 
total decay width B meson, r 5 (= l/r^), we use the world average B life-time values 
(combined results from ALEPH, CDF, DELPHI, L3, OPAL and SLD) |2I]: 

Tpo = 1.546 ±0.021 ps , 

rs+ = 1.647 ±0.021 ps . (44) 

To compare theoretical results with experimental data, as well as to determine con¬ 
straints on the effective number of colours, ^ |;] 4 g form factors and the CKM ma¬ 
trix parameters, we shall use the experimental branching ratios collected by CLEO 
BELLE jHni 1^ IS2] and BABAR jSSl IMj factories. All the experimental values are sum¬ 
marized in Table m 

In order to determine the range of ^ which is allowed by experimental data, we have 
calculated the branching ratios for B^ —>■ p'^tt^, —>■ p^7r°, R° —>■ p^tt"^, and R° ^ p°7r°. 
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All the results are shown in Figs. El El El and |H1 for the corresponding branching ratios listed 
above. Results are plotted for models (1) and (2), since they involve different form factor 
values and thus show their dependence on form factors. As experimental data, we shall 
use three sets of data from the CLEO, BABAR and BELLE Collaborations, respectively. 
Since experimental branching ratios from CLEO are the most accurate, we shall use them 
to extract the range of . The other two, the BABAR and BELLE data, will give 
us an idea of the magnitude of the experimental uncertainties. It is clear that numerical 
results are very sensitive to uncertainties coming from the experimental data. Thus, the 
determination of the allowed range of will be done by using all the branching ratio 
results. 

Let us start with the decay processes B~ p^7i~ and B~ —>• p“7r°. In both cases, 
there is a large range of acceptable values for a^d the CKM matrix elements over 
which the theoretical results are consistent with experimental data from CLEO, BABAR 
and BELLE. For B~ the lack of data does not allow us to determine the range. 

However, experiment and theory are consistent in both cases. For B~ P^t^~■, the models 

show considerable variation even though they are all consistent with the experimental 
data. Numerical results for models (1, 3) and (5) are close, so are those for models (2) 
and (4). We emphasise that the effect of p —a; mixing on the branching ratio B^ 
can be as large as 30%. As regards B^ p~7r~^ and —>• p°7r°, the results and conclusions 
are different from those for B~ —> p^7r~. If we look at the branching ratio for B^ —>■ p^7^^, 
only models (2) and (4) are consistent with experimental data over a large range of ^ 
whereas models (1,3) and (5) are not. The strong sensitivity to the results in that case 
comes from the fact that the decay branching ratios for B^ —> p^ir^ depend on form 
factors more sensitively, because in this case only one form factor, Fi(fc^), is involved. In 
all the other cases, the amplitudes depend on both Fi{k^) and Ao^k^). Therefore these 
branching ratios are less sensitive to the magnitude of the form factors. Finally, for the 
branching ratio BR{B^ —> utn^) plotted in Fig. El all models give theoretical results 
consistent with experimental data. Once again, the difference observed between models 
(1) and (2) mainly comes from the form factor (i.e. from the pion wave function 

used). Our analysis shows that models (1,3) and (5) cannot give results consistent with 
all experiments and have to be excluded. 

To remove systematic uncertainties coming from experimental results, one can cal¬ 
culate the ratio between two branching ratios for B decays. In the present case (with 
the data available), the ratio, is between BR{B^ p^n^) and BR{B^ p^7T^). 
Results are shown in Fig.^l We observe that the ratios differ totally from each other for 
models (1,3) and (5) and models (2) and (4). Since models (1,3) and (5) have already 
been excluded, we will use models (2) and (4) for the determination of the range for . 
If we just include tree contributions in the decay amplitudes, becomes independent of 
the CKM matrix elements. Penguin contributions lead to a relatively weak dependence 
of Rt^ on the CKM matrix elements. By comparing numerical results and experimental 
data, we are now able to extract a range for -which is consistent with all the results. 
To determine the best range of ^ -we select the values of -which are allowed by all 
constraints for each model. Finally, after excluding models (1,3) and (5) for the obvious 
reasons mentioned before, we can now fix the upper and the lower limit of the range of 
Nt" (TableEl). We find that should be in the range 1.09(1.11) < < 1.68(1.80) 
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for = 0.3(0.5). Comparing with our previous study, the current range of is 

consistent but smaller than the previous one. 


6 Summary and discussion 

The hrst aim of the present work was to compare theoretical branching ratios for 

p^7r°, B^ —>■ and B^ p°7r° with experimental data from the CLEO, 

BABAR and BELLE Collaborations. The second was to apply recent values of the CKM 
matrix elements, e.g. A,X,p and p, to study direct CP violation for B decay such as 
B^ — >■ p^{uj)ti^ — >■ 7 r’''7r“7r‘^, where the p — uj mixing mechanism must be included. The 
advantage of including p — uj mixing is that the strong phase difference which is necessary 
for direct CP violation, is large and rapidly varying near the uj resonance. As a result, 
the CP violating asymmetry, acp, reaches a maximum, amax, when the invariant mass of 
the TT+TT” pair is in the vicinity of the uo resonance and sin 5 = +1 at this point. 

In our approach, we started from the weak effective Hamiltonian where short distance 
and long distance physics are separated and treated by a perturbative approach (Wilson 
coefficients) and a non-perturbative approach (operator product expansion), respectively. 
One of the main uncertainties introduced in our calculation comes from the hadronic 
matrix elements for both tree and penguin operators. We treated them by applying a naive 
factorization approximation, where N^f is taken as an effective parameter. Although this 
is clearly an approximation, it has been pointed out |2Sj that it may be quite reliable in 
energetic weak decays such as B ^ pn. 

We have investigated the direct CP violating asymmetry in the B decay: B^ 
7T~^7r~7T^. We found that the CP violation parameter, acp, is very sensitive to the param¬ 
eters p and r] in the CKM matrix, and also to the magnitude of the form factors appearing 
in the hve phenomenological models we investigated. We have calculated the maximum 
asymmetry, amax, as a function of the effective parameter, , with the limiting values 
of the CKM matrix elements. We found that the CP violating asymmetry, Qmax, can 
vary from —37% to —84% over all the models (1,2, 3,4, 5). As we already suggested in 
a previous study I2HI. the ratio between the asymmetries for limiting values of the CKM 
matrix elements is mainly governed by p. Previously, we found a ratio equal to 1.64 where 
the CKM values used were the following: A = 0.815, A = 0.2205, 0.09 < p < 0.254, and 
0.323 < p < 0.442. In the present work, we found for the same decay, a ratio equal to 
1.30. The more accurate value for p has reduced uncertainties on both the CP violating 
asymmetry and the ratio, r(i?^ — )/r(B» ^ p» 7r°). 

Moreover, we stressed that without the p — uj mixing mechanism, the CP violating 
asymmetry, acp (which is proportional to both sin5 and r), is small since in that case 
either sin 5 or r is small. In the allowed range of , we also found that the sign of 
sin^ is always positive. Therefore, by measuring acp, we can remove the phase mod(7r) 
ambiguity which occurs in the usual method for the determination of the CKM unitarity 
angle a. 

We have calculated branching ratios for B^ ^ B^ —> B^ —>■ and 

BO pO^o compared the results with experimental data coming from the CLEO, 
BABAR and BELLE Collaborations. We have shown that for models (2) and (4) there 
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is a range for 1.09(1.11) < < 1.68(1.80), in which theoretical resnlts are 

consistent with experimental data. Models (1, 3) and (5) are exclnded since the form factor 
Fi(fc^) in these models cannot prodnce resnlts consistent with experiment. For a deeper 
investigation into this problem, some resonant and non-resonant contribntions pni EH 
which may carry bigger effects than expected in the calcnlation of branching ratios in pvr 
may have to be considered serionsly. 

With more accnrate CKM matrix elements valnes, e.g. p and p, we are able to give 
more precise CP violating asymmetries, and the main nncertainties remaining are from 
the factorization [H] approach and the hadronic decay form factors. In the fntnre one may 
hope to nse QCD factorization to replace the effective parameter, , and hence to 
provide a more reliable treatment of non-factorizable effects. With regard to form factors, 
we have shown that some models for the 5 —>■ vr transition are not consistent with the 
experimental branching ratios. We expect that onr predictions will provide nsefnl guidance 
for future investigations in B decays. We look forward to even more accurate experimental 
data from our experimental colleagues in order to further constrain our theoretical results 
and hence, to further advance the determination of the CKM parameters p and p and our 
understanding of CP violation within or beyond the Standard Model. 
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Figure 1: CP violating asymmetry, acp, for for k‘^/ml = 0.3, N^f = 

1.09(1.68) and limiting values of the CKM matrix elements for model (1): solid line (dotted 
line) for = 1.09 and max(min) CKM matrix elements. Dashed line (dot-dashed line) 
for N^f = 1.68 and max(min) CKM matrix elements. 



Figure 2: CP violating asymmetry, acp, for —>■ 7r’''7r“7r°, for k‘^/ml = 0.5, N^f = 
1.11(1.80) and limiting values of the CKM matrix elements for model (1): solid line (dotted 
line) for N^f = 1.11 and max(min) CKM matrix elements. Dashed line (dot-dashed line) 
for = 1.80 and max(min) CKM matrix elements. 
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Figure 3: sin5 as a function of , for —>■ for = 0.3(0.5) and 

for model (1). The solid (dotted) line at sin5 = +1 corresponds the case IIp^ = 
(—3500;—300), where p — oj mixing is included. The dot-dashed (dot-dot-dashed) line 
corresponds to IIp^ = (0; 0), where p — oj mixing is not included. 



Figure 4: The ratio of penguin to tree amplitudes, r, as a function of , for —>■ 
7i~^7r~7T^, for = 0.3(0.5), for limiting values of the CKM matrix elements [piP] 

max(min), for Ilp^ = (—3500; —300)(0, 0), (i.e. with(without) p — oj mixing) and for 
model (1). Figure 4a (left): for = (0;0), solid line (dotted line) for k‘^/m\ = 0.3 
and {piP) max(min). Dot-dashed line (dot-dot-dashed line) for k‘^/m\ = 0.5 and {p^p) 
max(min). Figure 4b (right): same caption but for Ilp^ = (—3500; —300). 
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Figure 5: Branching ratio for ^ for models 1(2), k‘^/m% = 0.3 and limiting 
values of the CKM matrix elements. Solid line (dotted line) for model (1) and max(min) 
CKM matrix elements. Dot-dashed line (dot-dot-dashed line) for model (2) and max(min) 
CKM matrix elements. Horizontal dotted lines: CLEO data; horizontal dashed lines: 
BABAR data; horizontal dot-dashed lines: BELLE data. 



Figure 6: Branching ratio for for models 1(2), k^/m% = 0.3 and limiting 

values of the CKM matrix elements. Solid line (dotted line) for model (1) and max(min) 
CKM matrix elements. Dot-dashed line (dot-dot-dashed line) for model (2) and max(min) 
CKM matrix elements. Same notation for experimental data as in Fig. El 
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Figure 7: Branching ratio for for models 1(2), k^/m% = 0.3 and limiting 

values of the CKM matrix elements. Solid line (dotted line) for model (1) and max(min) 
CKM matrix elements. Dot-dashed line (dot-dot-dashed line) for model (2) and max(min) 
CKM matrix elements. Same notation for experimental data as in Fig. El 



Figure 8: Branching ratio for —> p°7r° for models 1(2), /m\ = 0.3 and limiting 

values of the CKM matrix elements. Solid line (dotted line) for model (1) and max(min) 
CKM matrix elements. Dot-dashed line (dot-dot-dashed line) for model (2) and max(min) 
CKM matrix elements. Same notation for experimental data as in Fig. El 
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Figure 9: Branching ratio for —>■ (uvr^ for models 1(2), /m\ = 0.3 and limiting 

values of the CKM matrix elements. Solid line (dotted line) for model (1) and max(min) 
CKM matrix elements. Dot-dashed line (dot-dot-dashed line) for model (2) and max(min) 
CKM matrix elements. Same notation for experimental data as in Fig. El 



Figure 10: The ratio of two pn branching ratios versus for models 1(2) and for 

limiting values of the CKM matrix elements: solid line (dotted line) for model (1) with 
max(min) CKM matrix elements. Dot-dashed line (dot-dot-dashed line) for model (2) 
with max(min) CKM matrix elements. Same notation for experimental data as in Fig. El 
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q‘^ /ml = 0.3 

q‘^/ml = 0.5 

c*; 

-0.3125 

-0.3125 


+1.1502 

+1.1502 

^3 

+2.433 X 10"^ + 1.543 x lO’^i 

+2.120 X 10-2 + 2.174 X lO-^i 


-5.808 X 10-2 - 4.628 x lO-^z 

-4.869 X 10-2 - 1.552 x 10-2z 

a 

+1.733 X 10-2 + 1.543 x lO-^i 

+1.420 X 10-2 + 5.174 X lO-^i 

C'e 

-6.668 X 10-2 - 4.628 x lO-^i 

-5.729 X 10-2 - 1.552 x 10-2^ 


-1.435 X 10-^ - 2.963 x lO-'^z 

-8.340 X 10-3 - 9.938 x lO-^z 

^8 

+3.839 X 10-^ 

+3.839 X 10-^ 


-1.023 X 10-2 - 2.963 x lO-'^z 

-1.017 X 10-2 - 9.938 X lO-'^z 

C'lo 

+1.959 X 10-3 

+1.959 X 10-3 


Table 1: Effective Wilson coefficients for the tree operators, electroweak and QCD penguin 
operators (see Refs mma). 



^Ao 

h 

mAo 

mi 

do{di) 

bo{bi) 

model (1) 

0.280 

0.290 

5.27 

5.32 



model (2) 

0.340 

0.625 

5.27 

5.32 



model (3) 

0.280 

0.290 

5.27 

5.32 



model (4) 

0.340 

0.625 

5.27 

5.32 



model (5) 

0.372 

0.305 



1.400(0.266) 

0.437(-0.752) 


Table 2: Form factor values for B ^ p and R —vr at = 0. 



n:il = 1.09(1.11) 

NgL = 1.68(1.80) 

model (1) 

Pmax’) Pmax 

-55(-41) 

-65(-51) 

Pmin ? Pmin 

-72(-55) 

-80(-65) 

model (2) 

pmax-) Pmax 

-63 (-48) 

-71(-56) 

pmin 1 Pmin 

-78(-62) 

-84(-69) 


model (3) 


pmax-) Pmax 

-56(-41) 

-65(-51) 

pmin’) Pmin 

-72(-55) 

-80(-69) 


model (4) 


Pmax •) Pmax 

-64(-48) 

-71(-57) 

pmin-) Pmin 

-79(-62) 

-84(-69) 

model (5) 

Pmax 5 Pmax 

-51(-38) 

-58(-44) 

pmin ’) Pmin 

-63(-51) 

-72(-60) 


Table 3: Maximum CP violating asymmetry amax{%) for for all models, 

limiting (upper and lower) values of the CKM matrix elements, and k‘^/ml = 0.3(0.5). 
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CLEO 

BABAR 

BELLE 


10.4+j|;^ ± 2.1" 

24±8±3^ (< 39)^ 

8.04:44" (< 28.8)1 

H- 

o 

<43^1 

— 

— 

p'^n^ 

27.6tf,t ± 4.2" 

28.9 ±5.4 ±4.3* 

20.84:44" (< 35.7)1 

P^TT^ 

1.6^^;^±0.8' (< 5.5)^ 

< 10.6^1 

< 5.31 

BR{p^7t^) 

BR(pOn±) 

2.65 ± 1.9 

1.20 ±0.79 

2.60 ± 1.31 

- ± — 

OOTT 

11.3^4.9 ± 1-4^ 

6.6^4 ±0.7" 

4.24 8 ± 0.5" 


Table 4: The branching ratios measured by CLEO, BABAR and BELLE factories for 
decays into pn in unit of 10“® (see Ref in text). Experimental data*, preliminary results 
fit* and upper limit^. 


R —> pvr 

<[ [ with mixing 

model (2) 

1.09;1.63(1.12;1, 

.77) 

model (4) 

1.10;1.68(1.11;1, 

.80) 

maximum range 

1.09;1.68(1.11;1, 

.80) 

minimum range 

1.10;1.63(1.12;1, 

.77) 


Table 5: Best range of determined for k‘^/ml 


0.3(0.5) and for all B ^ pn decays. 
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